Strong correlations and the anisotropy of acceptor states in insulating La2-xSra;Cu04 
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We use the Green's function formalism to discuss the role of strong correlations to the spatial 
structure of acceptor states doped into a two-dimensional Mott-Hubbard antiferromagnetic insulator. 
When the scattering between doped carriers, at the nesting wave vector Q — {n, vr), is strong enough 
to produce a momentum dependent scattering rate, Fk, the corresponding acceptor states become 
spatially anisotropic. As an example, we calculate the spatial structure of an acceptor state bound 
to an attractive two-dimensional Dirac delta potential, for a simple form of Fk- We then discuss 
the role of such spatial anisotropy for the understanding of an apparent discrepancy between low 
temperature transport data and photoemission spectra in lightly doped La2-a:Sr2:Cu04. 

PACS numbers: 78.30.-j, 74.72. Dn, 63.20.Ry, 63.20.dk 



The quantitative description of hopping conductance 
in extrinsic semiconductors relies on what is called the 
effective mass approximation^. Once the band structure 
of the host material is known, the effective masses on 
top of the valence band, for acceptors, or bottom of the 
conduction band, for donors, can be calculated, and a 
Schrodinger equation for the trap potential, generated 
by the impurity ions on the doped carriers, is then writ- 
ten down^. The bound states associated with the nega- 
tive energy solutions to this equation correspond to local- 
ized states and, usually, the envelope function at the low- 
est energy state is spatially isotropic, symmetric, parity 
even, just like the Hydrogen atom Is state. The Bohr ra- 
dius of the bound state (or inverse effective mass) controls 
the exponential decay of the envelope wave function, the 
overlap between two spatially separated impurity states, 
and directly affects the hopping conductance. 

Undoped cuprates, such as La2Cu04, are far from be- 
ing classified as extrinsic semiconductors. These insu- 
lating materials actually belong to the class of strongly 
correlated electron systems, and exhibit long ranged an- 
tiferromagnetic order in the ground state. Nevertheless, 
it has been found that, when doped with a small amount 
of carriers, La2_2:Sr2:Cu04, for example, exhibits metal- 
lic behavior at high temperature (albeit anomalous) al- 
ready for 1% of carriers^. Furthermore, La2_xSr2:Cu04 
also exhibits hopping conductance at low temperature, in 
a very similar fashion as described above. In particular, 
at X = 3% and x = 4%, both AC^ and DC^ data are 
consistent with hopping transport and point towards the 
existence of localized acceptor states. 

Contrary to the above tendency, recent angular re- 
solved photoemission spectroscopy (ARPES) experi- 
ments in La2_xSra;Cu04 revealed a certain degree of 
"metalicity" at the same doping regime discussed above^, 
and have, as a result, been interpreted as evidence for 
delocalization of the charge carriers. The observation of 
"quasiparticle" peaks in the form of small Fermi arcs for 
3% doped La2_a;Sr:ECu04 at T = lOK—, which according 
to transport should be insulating, resembles very much 
the full pocket structure observed in other cuprates at 
higher doping^, where the sample is (strangely) metallic. 
The conclusion was that, since ARPES shows "peaks" in 



reciprocal (momentum) space, the associated wave func- 
tion should be extended in real (coordinate) space, and 
not localized as suggested by hopping transport. 

As an attempt to address this apparent discrepancy, 
A. S. Alexandrov and K. Reynolds® suggested that 
such Fermi arcs observed with ARPES in insulating 
La2_2:Sr2:Cu04 could arise naturally from matrix ele- 
ments when the ejected electron leaves behind a hole at 
valence band tails, formed by the hybridization of valence 
band and impurity states^. The reasoning is actually rel- 
atively simple. While long-lived quasiparticles in a Fermi 
liquid have very well defined momentum states and thus 
give rise to sharp peaks in the ARPES spectrum, local- 
ized wave functions are written as a linear superposition 
of many different wave vectors, and, as a consequence, do 
not possess a well defined momentum state. For this rea- 
son, no peaks are to be expected in the ARPES spectrum 
corresponding to a bound state. For a sufficiently large 
envelope wave function, however, the ARPES matrix el- 
ements in reciprocal space will depend on which of the 
linearly combined wave vectors weigh a heavier contribu- 
tion to the spectral function. For the case of cuprates, 
for example, it is argued, in Ref. 's', that the wave vec- 
tors close to (77/2, tt/ 2) contribute a much higher spectral 
weight, which then falls rapidly away from this point^. 
For this reason, bound electrons ejected from band tails 
are able to produce Fermi arcs in the ARPES spectrum. 

A necessary condition for the Alexandrov and 
Reynolds idea to work is that the localized state should 
be anisotropic in real space?, being more elongated along 
the direction perpendicular to the Brillouin zone (B.Z.) 
faces, and being shorter along the direction parallel to 
i1^. In this work we extend the traditional effective mass 
approximation to include strong correlations, and we cal- 
culate the full spatial structure of acceptor states us- 
ing a phenomenological model to mimic light doping in 
La2^2;Sr2;Cu04. We argue that due to the presence of a 
nesting AF wave vector Q = (tt, tt) (in units of the lat- 
tice spacing a — 1) the localized states become squeezed 
exactly along the B.Z. faces and elongated along the so 
called nodal directions (along the direction perpendicu- 
lar to the B.Z. faces), providing the missing ingredient to 
complete the Alexandrov and Reynolds analysis. 
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We begin by reviewing the effective mass approxima- 
tion for the calculation of acceptor/donor states in ex- 
trinsic semiconductors^. The effective mass tensor, m*^, 
is defined as (ti*^)"^ — 9^e(k)/9fcQ,fc^, where £(k) is the 
quasiparticle dispersion close to the bottom of the con- 
duction band (for donors) or top of the valence band (for 
acceptors) located at the i-th pocket. In case the effective 
mass tensor is diagonal, but with different components, 
m|J and m^, the time- independent Schrodinger equation 

for the localized wave function ip^{r) reads 



QuasiparticI; dispcnkm — no nestiQg 



Quasiparticle Dispmkm ^ nesting 
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where Vj. = d/dz, Vii 



+ V{r) + U{r) i/.*(r) V''(r), 

(1) 

{d/dx,d/dy), V{r) is the pe- 
riodic potential provided by the lattice, and U (r) is the 
trap potential felt by the doped carriers and provided by 
an impurity at the origin of the coordinate system. For 
large enough bound states, when the localization length 
is of the order of several lattice spacings, one can write 



(2) 



where _F*(r) ^ er"^!^ is usually an exponentially decay- 
ing envelope wave function with localization length 
at the i-th pocket and, according to Bloch's theorem. 



(/)*(k',r) — e ■''wi(;i(r), where u-^i(v) is a periodic func- 
tion with minimum at k' . Following Kohn and Luttinger 
(KL)^, each one of the [i = \ . . .N degenerate (from dif- 
ferent pockets) localized states are generally written as 



N N 
i=l 4=1 



E a^i^^ (r)0^(k\r), (3) 



where are A^-dimensional vectors determined by point 
group symmetry carrying the information about the sym- 
metries of the wave functions. We define 7 — m*j_/m'^^, 

and F^{r) ^ e~^^^' can be classified as isotropic 3-D, 
7 = 1, in which case r = X'^ + + Z'^ , or extremely 
adiabatic 2-D, 7 — > 0, in which case r = y/x^ +2/^- 

It has long been acknowledged that, to a very good 
extent, the physics of lightly doped cuprates is well cap- 
tured by the t ~ t' — t" ~ J modeU^, which includes 
not only nearest neighbor hopping, i, but also second 
neighbor hoppings, t' and t", for holes doped into a two- 
dimensional antiferromagnet with superexchange J. At a 
very small concentration, the Fermi surface for the added 
holes correspond to four, nearly two-dimensional hole 
pockets centered at the wave vectors ko = (±7r/2, ±7r/2) 
in the magnetic B.Z., as can be seen in Fig. [T}i). This 
is consistent with band structure calculations, as well as 
ARPES data^, and the two-dimensionality of the Fermi 
surface results from the huge effective along the direc- 
tion perpendicular to the Cu02 planes (or parallel to 
the c-axis), ^ mj|, where TOjj is the planar effective 
mass. The above large effective mass anisotropy drives 



FIG. 1: (Color online) Evolution of the quasiparticle spec- 
trum: a) (left) fully coherent hole-like Fermi pockets^, lo- 
cated at ko = (±7r/2, ±7r/2) in the magnetic B. Z.; b) (right) 
additional scattering at the nesting antiferromagnetic wave 
vector^'^, Q = {-k^-k), and the hot spots along the magnetic 
B.Z. faces (effect of a k-dependent scattering rate Pk). 



the system into the extreme adiabatic limit discussed 
above, 7 — > 0. Self consistent Born approximation for 
the t — t' — t" — J model has also determined that, at 
the very top of the valence band, the pockets are nearly 
circular, with equal effective masses along the magnetic 
B.Z. faces and diagonals, or orthorhombic (a, h) axis, and 



thus we can set m* 



In this case 



F[v) 



(4) 



where r — X'^ + y-^ and ^ ~ Such large effective 

mass anisotropy, m*^ ^ mj| is consistent with the large 

resistivity anisotropy in these materials^, Pc/ Pab ^ 10^. 

However, a 2D isotropic wave function is not consis- 
tent with neither the ARPES response for very lightly 
doped La2_a;Sra;Cu04 nor with the AC and DC trans- 
port measurements at low temperature, in the hopping 
regime, where the knowledge of the precise shape of 
the localized state turns out to be crucial. These ex- 
periments have revealed that, besides the perpendicular 
anisotropy discussed above, m^/mj| ^ I, an additional 
in-plane anisotropy is necessary to account for the differ- 
ent infrared absorption spectra observed along the two 
orthorhombic directions in La2--a;Sr2;Cu04^. In what fol- 
lows, we shall argue that, even more important than the 
anisotropy caused by spiral correlations considered by V. 
Kotov and O. Sushkov in Ref. [ll|, is the anisotropy com- 
ing from strong correlations, in the form of a momentum 
dependent scattering rat o^^i^^ . This will force us to ex- 
tend the usual effective mass approximation, used in Eq. 
P^, in order to incorporate such nontrivial effect. 

On general grounds, we shall assume that the origi- 
nal quasiparticle dispersion of the noninteracting, J = 0, 
problem, the t~t' ~t" model, have the pocket like struc- 
ture at ko = (±7r/2, ±7r/2) shown in Fig. [TJi). For J 7^ 0, 
in turn, short wavelength antiferromagnetic fluctuations, 
at low temperature and near half-filling, cause the en- 
hancement of the quasiparticle scattering at the ordering 
wave vectorJi, see Fig. [T}d). The presence of such hot 
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spots produces, in turn, a momentum dependent scatter- 
ing rate (or inverse quasiparticle lifetime}!^, which pro- 
vides us with an important source of anisotropy. 

For siniphcity. let us consider the problem of a hole 
under the influence of an attractive 2D-delta potential^. 
According to the usual effective mass approximation, the 
envelope wave function F{r) can be obtained from 



' 2mr. 



U{v) 



F{r) = E F(r) 



where 



t/(r) 



-9- 



^(r-ro) 



(5) 



(6) 



is an attractive (5-potential that traps, with strength g, a 
hole around an impurity ion at the origin. We shall be 
interested in the solutions to this equation for r larger 
than the cutoff Tq. Using the Fourier transform 



(7) 



the equation becomes 



(2^ 



2mt 



-E 



(8) 

We now multiply the whole equation by (27r)^ e"^ '', and 
integrate over r to obtain 



/(q) 



47r^ffo F(ro) 
e(q) + SB 



= 90 



fjro) 
e(q) + £B ' 



where we defined = g/r^^, /(''o) — 47r^F(ro), and 
£(q) = /i^q^/2m|j, is the dispersion in the effective mass 
approximation. Here we have used, as the negative en- 
ergy (bound state) solution, —E = eb — fi^K^ /2m*^, with 
K = 1/^ playing the role of inverse localization length ^. 
We recognize Eq. © as a particular case of the more 
generic Bethe-Salpeter equation for a potential U (r) 



/(q) 



47r2 



e(q) +eB J (27r) 



U{k' - q) /(k'), (10) 



with C/(k) being, as usual, the Fourier transform of U{r). 

The result presented in ([9]) is the exact solution to the 
equation ^ in which all wave vector are treated on equal 
footing. In particular, the term fi,^V'^/2mJ| is isotropic 
and results from the nearly perfect circular shape of the 
hole pockets in lightly doped La2-a;Sra;Cu04. For this 
reason, /(q) depends only on q^ and the resulting enve- 
lope function F{r) is also isotropic. In order to introduce 
the momentum dependence of the scattering ratei^, Ffc, 
a quantity that is related to the quasiparticles spectral 
function and which does not appear in ([SJ , we will make 
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FIG. 2: (Color online) Anisotropic scattering rate (red 
solid line) as a function i^. Inset: is the angle between 
the wave vector q (black solid arrow) and the magnetic B.Z. 
face (black dashed line). At the hot spots, = 0, tt, the 
scattering is maximal, while at the nodal direction, i\> = 7r/2, 
it is mininial^^. 



use of the Green's functions formalism, and its spectral 
representation. We perform a Stieltjes transform to write 



G(q,L. = 0) 



1 



e(q) + SB 



dE^^, (11) 



where the spectral function A(q, E) is, for such strongly 
correlated electron system, given by 



Aiq,E) = 



TT (iJ - e(q))2 + r2 



+ A,„c{q,E), (12) 



(9) where Zq is the quasiparticle weight, Fq is the momen- 
tum dependent scattering rate, and Ainc is the incoherent 
(or multiparticle) part of the spectrum. For a weakly in- 
teracting system, a Fermi liquid, where Fq = Fq ^> 0, 
one ends up with A{q, E) = ZqS{E — e(q)) -I- Ai„c(q, E). 

Since the nesting is strongest for the points where the 
pockets intersects the B.Z. faces, we can, without loss of 
generality, consider the phenomenological formula 



Fq |cOs((/>) 



(13) 



where </> is the angle between the wave vector q, centered 
at a given pocket, and the B.Z. face. This simple form of 
the scattering rate produces larger widths for the spectral 
lines along the B.Z. faces, the hot spots at = 0, tt, very 
small widths along the nodal directions, (j> = Tr/2 (see 
Fig. [2]), and is consistent with the results from Ref. \13i 

For calculating bound states, it is enough to keep 
only the coherent part of the spectrum, defined as 
Acoh {<i, E) = A{q, E) — Ai„c(q, E), and thus we write 



F(r) = 5o/(0) 



(2^)2 



dE l^ AcoiM^E). (14) 



E + eb 



We now make the following approximations. We use ex- 
plicitly that Zq = Zq, that is, the quasiparticle weight is 
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weakly momentum dependent at low doping, and we use 
explicitly the angular dependence of the scattering rate, 
Tq = r^. We evaluate this integral by closing a contour 
in the upper half plane and use the residue theorem 



F{r) = 9of{ro)Zo 



(2^)2 



e(q) + SB 



(e(q)+£B)2 



p2 ■ 

^0 



(15) 

where we have used explicitly that Zq — Zq and Fq = F^. 
Notice that in the limit where F^ — > and Zq — > 1 the 
above equation reduces to the result obtained in Q. 

Before performing the angular integration, we must 
recall that, while cj) is the angle between q and the B.Z. 



face, the angle that appears in q • r = qr cos 9, is the 
angle between q and r. So, if we define as (p the angle 
between r and the B.Z. face at a given pocket, we end 
up with = (j) — tp. The complete solution to the bound 
state problem can be finally presented. We approximate 
F^ < es, and write F^{r) = Fo{r) + SF^{r), where 



Fo(K,r) = gof{ro)Zo^Ko{Kr), 



(16) 



is the unperturbed (Fq = 0) isotropic (iy9-independent) 
contribution to the wave function, while the perturbation 
(and source of anisotropy in (p) becomes 



J 



2(TOp2 j 2 (2 - K'^r'^KoiKr)) cos (2(^) - KrKi{Kr) (^^r^ + (4 + kV^) 



cos 



(2^)) I 



(17) 



r 



y o\ 



-10 1! 



Density probability at (;r/2,-7r/2) 



Density probability at (jr/2 jr/2) 




y 



-10 -5 5 




-10 -5 5 10 



FIG. 3: (Color online) Density plot of the amplitude probabil- 
ity for the envelope wave function F(r), at the two nonequiva- 
lent pockets: a) (left) (7r/2, -7r/2); and b) (right) {n/2,n/2). 
The axis a and b are the orthorhombic axis. 



where Kq and Ki are modified Bessel functions of the 
second kind and, from its definition, (p is given by 
p = arctan[(a; =F y)/{x ± y)], for the pockets centered 
at (7r/2, ±7r/2), see Fig. [31 We see that, as required by 



Alexandrov and Reynolds^, the spatial structure of the 
acceptor state is such that the envelope wave function is 
elongated exactly along the nodal directions, or perpen- 
dicular to the B.Z. faces, and squeezed along the direc- 
tions parallel to the B.Z. faces. 

We have extended the traditional effective mass 
approximation^i to include the effects from strong cor- 
relations, reflected in the anisotropy of the quasiparticle 
scattering rate^^, Fk (or inverse quasiparticle lifetime). 
Using the Green's function method, and its spectral rep- 
resentation, we were able to calculate the full structure 
of an acceptor state trapped by a simple two-dimensional 
(5-potential, used to mimic the problem of lightly dop- 
ing a Mott-Hubbard antiferromagnetic insulator. Our 
results provide the missing ingredient to the Alexandrov 
and Reynolds analysis^, which bridges the gap between 
low temperature transport^ and ARPES^ data. 

The author acknowledges discussions with A. Alexan- 
drov, G. Blumberg, R. Gapaz, A. Ghernyshev, B. Koiller, 
and O. Sushkov. 
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